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Abstract 

The risk of internal exposure to ^^^Cs, ^^"^Cs, and ^^^1 is of great public concern after the accident at the Fukushima-Daiichi 
nuclear power plant. The relative biological effectiveness (RBE, defined herein as effectiveness of internal exposure relative 
to the external exposure to y-rays) is occasionally believed to be much greater than unity due to insufficient discussions on 
the difference of their microdosimetric profiles. We therefore performed a Monte Carlo particle transport simulation in 
ideally aligned cell systems to calculate the probability densities of absorbed doses in subcellular and intranuclear scales for 
internal exposures to electrons emitted from ^^^Cs, ^^"^Cs, and ^^^1, as well as the external exposure to 662 keV photons. The 
RBE due to the inhomogeneous radioactive isotope (Rl) distribution in subcellular structures and the high ionization density 
around the particle trajectories was then derived from the calculated microdosimetric probability density. The RBE for the 
bystander effect was also estimated from the probability density, considering its non-linear dose response. The RBE due to 
the high ionization density and that for the bystander effect were very close to 1, because the microdosimetric probability 
densities were nearly identical between the internal exposures and the external exposure from the 662 keV photons. On the 
other hand, the RBE due to the Rl inhomogeneity largely depended on the intranuclear Rl concentration and cell size, but 
their maximum possible RBE was only 1.04 even under conservative assumptions. Thus, it can be concluded from the 
microdosimetric viewpoint that the risk from internal exposures to ^^^Cs, ^^"^Cs, and ^^^1 should be nearly equivalent to that 
of external exposure to y-rays at the same absorbed dose level, as suggested in the current recommendations of the 
International Commission on Radiological Protection. 
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Introduction 

The risk of internal radiation exposure is of great public concern 
after the accident at the Fukushima-Daiichi nuclear power plant 
[1,2]. This is partially because the risks from internal exposure can 
differ from those from external exposure to y-rays at the same 
absorbed dose level, i.e., their relative biological effectiveness 
(RBE, defined herein as effectiveness of internal exposure relative 
to the external exposure to y-rays) is not always 1. For example, 
there is evidence that RBEs for the intake of cx, low-energy (3, and 
Auger-electron emitters are greater than 1 [3,4]. In general, the 
high ionization density around the trajectories of a particles and 
low-energy electrons as well as the inhomogeneous radioactive 
isotope (Rl) distribution in subcellular structures are considered to 
explain the higher RBE values. These are referred to hereafter as 
the track-structure and Rl-inhomogeneity effects, respectively. 

A number of studies have been carried out to estimate the RBE 
for the intake of oc, low-energy P, and Auger-electron emitters [5- 
7], as these results are well summarized in Report of the 
Committee Examining Radiation Risks of Internal Emitters [4]. 
In contrast, the RBE for the intake of ^^^Cs, ^^^Cs, and ^^^I (major 
contributors to the internal exposure dose from the nuclear 



accident in Fukushima) was not extensively discussed. This is 
because these RIs emit relatively high energy electrons and 
photons, and because the scientific community considers that their 
RBE is 1. Nevertheless, the public occasionally believes that the 
risk from internal exposure is much greater than that from 

137 1 34 131 

external exposure even for the intake of ^"^^Cs, ^""Cs, and ^"T, 
albeit no supportive scientific evidence. Such belief comes, at least 
in part, from the lack of a detailed analysis of the contribution of 
the track-structure and Rl-inhomogeneity effects to the RBE for 
the intake of these RIs, except for the Rl-inhomogeneity effect of 
^^^I [8]. The contribution of these effects may not be negligible, if 
these RIs are selectively located inside cell nucleus. An animal 
study has suggested that 12-21% of ^^^Cs are localized in cell 
nuclei [9]. 

We therefore set out to quantitatively analyze the contribution 
of the track-structure and Rl-inhomogeneity effects for the intake 
of ^^^Cs, ^"^^Cs, and ^^^I, and for this, a microdosimetric 
simulation was performed using the Particle and Heavy Ion 
Transport code System (PHITS) version 2.64 [10]. The RBE for 
internal exposure to these RIs was then derived as a function of 
cell size and the fraction of the Rl distributed in cell nuclei. In 
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addition, the RBE for the bystander effect (biological effect caused 
by signaling from irradiated to non-irradiated cells) [11,12] was 
also discussed based on the inhomogeneity of absorbed dose 
among cell nuclei, since the bystander effect may play an 
important role in the risk estimation of low-dose internal exposure 
owing to its non-linear dose response [13,14]. The results of the 
simulation, together with the maximum possible RBE from the 
dosimetric viewpoint, are presented in this paper. 

Materials and Methods 

Monte Carlo Particle-transport Simulation in Cell Systems 

We performed Monte Carlo particle transport simulations in 
ideally aligned cell systems that were internally exposed to 
electrons emitted from ^^^Cs, ^^^Cs, and ^"^^I, using the PHITS 
code. Electrons emitted from ^^^™Ba, which is a daughter isotope 
of ^^^Cs with a half-life of 2.552 min, were also considered in the 
simulation of the intake of ^^^Cs. Similar simulations were also 
performed for internal exposures to electrons from ^H, ot particles 
from ^^^Pu, and 662 keV mono-energetic photons, which is the 
dominant y-rays from ^^^Cs (in strict sense, 662 keV photons are 
emitted from the decay of ^^^"^Ba). The last simulation condition 
also represented the external exposure to y-rays, because the 
source location is not an important factor for photon exposure 
from the microdosimetric viewpoint. Thus, this condition served as 
the reference condition in this study, i.e., its RBE is equal to 1. The 
energy spectra of particles emitted from the RIs were taken from 
the International Commission of Radiological Protection (ICRP) 
Publication 107 [15]. 

The geometry of the simulations is shown in Fig. 1. Cells and 
their nuclei were assumed to be concentric spheres comprised of 
1 g/ cm^ of liquid water. They were placed in an 11x11x11 lattice 
structure, yielding 1,331 cells in the system. Cell nuclei were 
categorized into six groups according to the distance from their 
center to the origin of the cell system, L (Fig. 1). The cell system 
was infinitely surrounded by liquid water. To analyze the cell size 
dependence, we prepared 15 different cell systems by uniformly 
changing the radii of cells and their nuclei, denoted tq and Tn, 
respectively. The was changed from 3 to 7 |Im in 1 |Im steps, 
whUe fc was set to 1.5, 2, or 3 times larger than r^. For each cell 
system and RI source, the simulations were carried out four times 
by changing the RI localizations, where RIs were uniformly 
distributed in the cell nucleus, cytoplasm, extracellular space, or 
entire region of the central lattice. 

In the PHITS simulations, all types of radiation were 
transported down to 1 keV, below which particles stop and 
deposit their entire energy at their location, except for positrons 
that cause pair production. This local approximation is adequate 
for our simulation because the ranges of 1 keV electrons and oc 
particles are negligibly short compared to cell size. The probability 
density (PD) of the specific energy <^ in a cell nucleus, ^{^), was 
estimated by calculating the energy deposited in each cell nucleus i 
per source emission. In addition, the dose PDs in each cell nucleus 
i as a function of the lineal energy j for a site diameter of 1 jim, 
dl^), were also derived from the PHITS simulation. The 
terminology for these microdosimetric quantities are defined in 
the International Commission on Radiation Units and Measure- 
ments (ICRU) Report 36 [16]. For calculating the unique 
microdosimetric function of PHITS [17] was employed, which can 
directly determine microdosimetric PDs using a mathematical 
equation developed on the basis of track-structure simulations 
[18]. 



Data Analysis 

The calculated PDs of <e in each cell nucleus, ^(4, were used to 
estimate the RBE due to the Rl-inhomogeneity and bystander 
effects. On the other hand, the calculated dose PDs as a function 
of J, were used to estimate the RBE due to the track- structure 
effect by combining the Q{y) relationship, which is the radiation 
quality factor expressed as a function of j for a site diameter of 
1 |Im defined in ICRU Report 40 [19]. For this estimate, we must 
convert the PDs for each cell nucleus to their mean value for all 
cell nuclei, including the contribution from cell nuclei located 
outside the lattice structure. The procedure for this conversion is 
given below. 

The mean number of cell nuclei categorized in group j having 
^>0, can be calculated by 



Jz>0 



(1) 



where / denotes the indicator function having a value of 1 only 
when —j, and gi is the group index for cell nucleus i. The single- 
event PDs of specific energy in a cell nucleus, /i gj{'^), were then 
determined by 
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,5" 



-im 



Gj 



0 



forz>0 , 
forz = 0. 



(2) 



Note that fi{z) was defined as the PD of z deposited in a single 
event in ICRU Report 36 [16], i.e., the contributions from non-hit 
targets are not included in fi{z)- The frequent mean specific energy 
in cell nuclei categorized in group j having ^>0, ^f,g/5 can be 
calculated by 



^F,G/-= |z/l,G/-(^)(iz. 



(4) 



To determine the contribution from cell nuclei located outside 
the lattice structure, which were categorized into group 7, we 
assumed that the single-event PD of z in those outer cell nuclei is 
the same as that for the outmost cell-nucleus group in the lattice 
structure, i.e., /i^g7('^) =/i,G6('^)- Consequently, zf,G7=%G6- The 
validity of this assumption will be discussed later in this paper. 
Under this assumption, the mean number of the outer cell nuclei 
having ^>0 per source emission, JVoy, can be calculated from the 
ratio of the total to the mean deposition energies in the outer cell 
nuclei, which is written as 



A^G7 = 



^F,G6^N 



^F,G6^L 



(5) 



where E^^^ is the total deposition energy outside the lattice 
structure, and and are the masses of a cell nucleus and a 
lattice, respectively. In this study, the value of £'out was determined 
by the PHITS simulation. The single-event PD of ^ in a cell 
nucleus averaged over the entire system, /i ^-^J^z), can be calculated 
by 
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Figure 1. Geometry of the cell system assumed in the PHITS simulation. 
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/l,ave(^) = 



E NgAgj{z) 

^ ave 



(6) 



where jV^ve is the mean number of the cell nuclei having z>0 per 

7 

source emission. Namely, A/'ave = E ^Gj- 

7=1 

Similarly, the dose PD as a function of j averaged over the 
entire system, ^aveW? also estimated from the dose PD for each 
cell nucleus i, djj), obtained from the PHITS simulation. Except 
for cell nuclei outside the lattice structure, the dose PD for cell 
nuclei categorized in group j, dcjiy), can be calculated by 



1331 

i=i 



(7) 



where Z/ is the mean specific energy in cell nucleus i per source 
emission. The value of Z/ can be determined by 



zMz)dz. 



(8) 



Forj= 7, the dose PD was again assumed to be identical to the 
outmost cell-nucleus group in the lattice structure, i.e., 
^GyW — ^ggW- Therefore, d^^Jj) can be calculated by 



E NGjZY,GjdGj{y) 

.7to 

7 

E NGjZv,Gj 
7=1 



(9) 



RBE Estimation 

RBE for the Rl-inhomogeneity effect. To estimate the 
RBE, we assumed that radiation effects are only initiated by the 
ionization inside a cell nucleus; although there have been several 
lines of evidence that targeted cytoplasmic irradiation can induce 
biological effects [20,21]. The reason for introducing this 
assumption is that our primary purpose is to estimate the 
maximum possible RBE for the internal exposure to ^^^Cs, 
^^^Cs, and ^^^I, which are expected to be higher under this 
assumption; when RIs are localized in cell nuclei, the higher cell- 
nucleus dose directly results in higher RBE values. 

Under this assumption, the RBE for the Rl-inhomogeneity 
effect can be defined as the ratio of the mean specific energy in a 
cell nucleus to that in a lattice. When RIs are uniformly distributed 
in each lattice and equilibrium between the incoming and 
outgoing particle energies is established, the mean specific energy 
in a cell nucleus and in a lattice, Zn and Zl, can be calculated by 



^N= ^/l,ave(^)dz 



(10) 



and 



(11) 
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respectively, where Ej^ is the total energy deposited inside a lattice, 
which corresponds to the mean source energy emitted from the 
RIs. Note that this RBE is equal to 1 for external exposure as well 
as internal exposure for the intake of RIs without any microscopic 
localization tendency, i.e., those uniformly distributed in all 
subcellular structures. 

It should also be mentioned that RIs are inhomogeneously 
distributed inside the human body not only on the microscopic 
scale of subcellular structures but also on the macroscopic scale of 
organs and tissues. However, the influence of the RI inhomoge- 
neity on the macroscopic scale has already been taken into account 
in estimating the effective dose for internal exposure by 
introducing biokinetic models as well as the tissue weighting 
factor [22]. Thus, the RBE due to macroscopic RI inhomogeneity 
is not discussed in this paper. 

RBE for the track-structure effect. In this study, the 
influence of the track- structure was represented by the mean 
quality factor based on the Q^) relationship [19], Q^ve? which can 
be calculated by 



4. Bystander signals uniformly propagate over a certain distance; 

5. All cells within the propagation distance can receive a 
bystander signal irrespective of whether they are directly 
irradiated or not; and 

6. The fraction of cells receiving a bystander signal from a single 
signal-emitting cell is constant. 

Except for items 3 and 5, these hypotheses are similar to those 
adopted in Ref. [28], and their adequacy was discussed therein. 

In the second and third hypotheses, the fraction of signal- 
emitting cells in the radiation field having the mean absorbed dose 
D, Psi-D), can be calculated by 



Ps{D)- 



|[l-5(z)]/av 



(14) 



(12) 



where [l-6*(4] represents the triggering probability of a cell having 
its nucleus specific energy z, and^ve('^5^) denotes the averaged PD 
of z in the radiation field. Similar to the analyses given in Ref 
[30], the values of^vefc^) can be determined by numerically 
solving the equation 



The RBE for the track-structure effect was then obtained from 
the ratio of Q^^^ for internal exposure to the reference condition, 
which was represented by exposure to 662 keV photons as above. 

It should be noted that the Q[L) relationship [23] is more 
frequently used than the Q^) relationship in RBE calculations for 
radiological protection purposes [3,24]. However, Q(Z) is not 
designed to express the energy dependence of the RBE for 
electrons except for very low energy to maintain the relation as 
much simplicity as possible. Instead, the use of Qj) allows RBE 
estimation considering the energy dependence, because the 
difference in RBE between y-ray and low-energy X-ray exposures 
is distinguishable in the relation. Thus, we here employed Q(j). 

RBE for the bystander effect. The bystander effect is 
considered attributable to the inhomogeneity of absorbed dose on 
microscopic scales. Thus, it should be in close connection with the 
PDs of z in cells or cell nuclei. Many studies have been devoted to 
develop models for quantitatively describing the bystander effect 
[25 29], but none of the existing models fully characterized the 
radiation fields using such microdosimetric PDs. We therefore 
developed an original model for describing the bystander effect 
based on the Fakir model [28], introducing the PDs of z in cell 
nuclei to express the dose inhomogeneity. 

The foUowings are the hypotheses adopted in our model: 

1. A cell is affected by bystander effects (e.g., those manifested as 
gene mutations, chromosomal aberrations, and cell killing) 
when receiving a bystander signal. 

2. Bystander signals are emitted from irradiated cells triggered 
with a probability depending on irradiation conditions, but 
their strength is independent of these conditions; 

3. The probability that a cell is not triggered after irradiated with 
its nucleus specific energy z, S{z), can be expressed by the 
linear-quadratic model in the same manner as its survival 
fraction, namely 



S{z)=Qxp{-az-l]z^) 



(13) 



where a and (3 are parameters that depend on radiation types 
and cell lines; 



,aveV^/ 



-m 



(15) 



k\ 



y^,ave(-^)? 



where P{X{D)\1i) represents the Poisson distribution with an 
expected value X{P) that is the mean number of events 
contributing to the dose, i.e., X{D)=D/z^. The function fk^aveiz) 
denotes the PD of z for cell nuclei irradiated by k events, which can 
be obtained from the convolution of those for cell nuclei irradiated 
by 1 and k—l events, yi ave('C) ^nd fk-i,sive{z), respectively, expressed 
as 



yA:,ave(-^) " 



3 (^1 )A - 1 ,ave — ^1 )dzi . 



(16) 



Using fi^aveiz) obtained from Eq. (6), the numerical values of 
Aave(4 were iteratively calculated up to when the conditions of k> 
X{D) and P{A{D);k)<OmOl were satisfied, e.g., k=l36 for 
X{D)=100. 

Let JVq be the number of cells within the distance in which the 
bystander signals can propagate, and rj be the fraction of cells 
receiving a bystander signal from a signal-emitting cell. The 
probability that a cell does not receive the entire signal emitted, 
Pa^, is then given by the sum of the binomial probabilities, which 
is written as 



Nq / at \ 

PAiD)= E ( „ j [^s(/))]"[l -Ps(/))]^«-"(l -»?)". (17) 



According to Ref [28], the calculation of this sum yields 
[l-f]Ps(D)f\ The probability that a cell receives a bystander 
signal, Pb{P^), which is referred to as the bystander probability in 
this paper, can be determined by 
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Pb(^) = 1 - PAiD) = 1 - [1 - r,PsiD)f« 



(18) 



It should be noted that this bystander probability does not 
represent the fraction of cells that actually respond to the 
bystander signal. For example, the actual fraction of cells 
inactivated after receiving a cell-killing bystander signal is 
supposed to be 10~20%, because the survival fraction of 
bystander cells is generally saturated around 80—90%, even for 
high-dose irradiation [31,32]. The actual fraction, however, is 
unlikely relevant to irradiation conditions and is hence not 
important in the estimation of the RBE. Thus, we defined the 
RBE for the bystander effect as the ratio of the absorbed dose that 
yields the same Pb{J^) for the internal exposures and reference 
condition. 

Results and Discussion 

Probability Density 

Figure 2 shows the calculated PDs of z for each cell-nucleus 
group j, zfi,Gj{z), for the exposure to electrons that are emitted 
from ^^^Cs localized inside the cell nucleus. It should be 
mentioned that the PD of z is generally depicted in the form of 
zf[z) on a semi-logarithmic graph, because their integrated 
probabilities, 7(<^)d<^, can be directly estimated from the graph by 
eye. These data are for the case of the smallest (r^ = 3 jim, 
rc = 4.5 jlm) and largest (^^=7 jlm, re = 21 jim) cells. It can be 
found from Fig. 2 that the PDs were nearly independent of the 
distance from the source to the cell nucleus except for j= 1, the 
central nucleus in which the RIs were localized. This tendency can 
be seen in the data for other irradiation conditions, thereby 
verifying the assumption that fi^Giiz) —f\,Ge,{z)- The peak observed 
in the PD forj= 1 for the smallest cell was due to the full-energy 
absorption of the Auger electrons around 4 keV. 

Figure 3 shows the averaged PDs multiplied by the mean 
number of irradiated cell nuclei per source emission, N^-vcZfi^z.-v(li)^ 
for exposure to electrons emitted from ^^^Cs localized in either the 
cell nucleus, cytoplasm, or extracellular space for the median cell 
size (rN = 5jIm, rc=10|lm). The PDs agreed well with one 
another except for the very low specific energy region, where the 
PDs become larger when ^^^Cs was localized in the cell nucleus. 
This is because low-energy (3-rays and Auger electrons can deposit 
their energy inside a cell nucleus only when they are generated 
inside or very close to the nucleus. The RBE for the RI 
inhomogeneity is attributed to this difference, as discussed later. 

Figure 4 shows the averaged PDs of z and Figure 5 shows dose 
PDs of J, zfi,sive{z) andj^/aveW, respectively, for exposures to various 
sources uniformly distributed inside a lattice. The PDs for high- 
energy P emitters, i.e., ^^^Cs, ^^""^Cs, and ^^^I, agreed well with one 
another, and were similar to those for 662 keV photons. 
Conversely, the PDs for a emitter ^^^Pu and low-energy (3 emitter 
were shifted to higher z or y regions because of their higher 
stopping powers, particularly for ^^^Pu. This verifies that the 
absorbed dose distributions in subcellular and intranuclear scales 
are inhomogeneous for the intake of d and low-energy (3 emitters, 
when compared to high-energy P emitters producing the same 
mean absorbed dose. The difference in zf[,a.ve{z) andj^/aveW results 
in the difference in the RBE for the bystander and track-structure 
effects, respectively, as discussed in the next section. 

RBE Estimation 

RBE for the Rl-inhomogeneity effect. Figure 6 shows the 
RBEs for the Rl-inhomogeneity effect as a function of radius of the 
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Figure 2. Dependence of the PDs of zon the cell-nucleus group 

/ These data represent the exposure to electrons emitted from ^^^Cs 
localized inside the cell nucleus. Panels A and B show the data for the 
smallest and largest cell sizes, respectively. 
doi:10.1371/journal.pone.0099831.g002 
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Figure 3. Dependence of the /Vave-weighted PDs of z on the RI 
localization tendency. Data represent the exposure to electrons 
emitted from ^^^Cs localized inside the cell nucleus, cytoplasm, and 
extracellular space, respectively, for a median cell size (rN = 5|im, 
rc = 10 |im). 

doi:10.1371/journal.pone.0099831.g003 
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Figure 4. Averaged PDs of z for exposure to various sources 
uniformly distributed inside a lattice. Data are for the median cell 
size (rN = 5 jim, rc = 1 0 |im). 
doi:l 0.1 371 /journal. pone.0099831 .g004 

cell nucleus. When the RIs is localized in cell nuclei, the RBE 
increases with increasing radii of cells and their nuclei due to the 
additional deposition energy from the low-energy particles emitted 
inside the cell nucleus. This tendency can be expected from the 
PDs shown in Fig. 3. The RBE was particularly high for and 
^^^Pu due to the shorter range of the emitted particles (low-energy 
(3 and CL particles, respectively). The RBE for ^^"^Cs was slightly 
higher than that for ^^^Cs and ^^^I because of its lower mean 
energy of the emitted P-rays. 

On the other hand, the RBE was less than or close to 1 when 
the RIs were localized in the cytoplasm or extracellular space 
because lower energy particles cannot reach the cell nucleus from 
such extranuclear compartments. Moreover, the RBE was nearly 
zero for localized in the extracellular space, thereby indicating 
that the (3-rays emitted from can deposit their energy inside the 
nucleus only when is incorporated into the cell. The RBE was 
occasionally greater than 1 for ^"^^Pu localized in the cytoplasm 
because some of the emitted a particles create a Bragg peak inside 
the cell nucleus. It should be noted that the RBE was always equal 
to 1 when RIs were uniformly distributed inside a lattice, although 
this was not shown in the figure. 

RBE for the track-structure effect. Table 1 summarizes 
the RBE for the track-structure effect estimated from ^aveW in 
combination with the Qj) relationship. These data represent the 
mean values and their standard deviations for all cell sizes and RI 
localizations analyzed in this study. Note that the standard 
deviations were mostly attributed to the differences in the 
calculation conditions because statistical uncertainties of the 
Monte Carlo particle transport simulations were negligibly small. 
Thus, the small standard deviations indicate that the RBE for the 
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Figure 5. Averaged dose PDs of / for exposure to various 
sources uniformly distributed inside a lattice. Data are for the 
median cell size (rN = 5 |im, rc=10 |im). 
doi:1 0.1 371/journal.pone.0099831 .g005 

track-structure effect is not sensitive to the size of cells or the 
localization of RIs. As expected from Fig. 5, the estimated RBE 

137 1 131 3 

was very close to 1 for Cs, Cs, and I, and those for H and 
^^^Pu were slightly and much larger than 1, respectively. 

It should be mentioned that the average quality factors, Q^^e? 
were smaller than the data given in Table 1 by a factor of 0.69, 
which is the value of Q^^ve for the reference condition (662 keV 
photon exposure). This is because the reference radiation of the 
Q(f) relationship was set to low-energy X-rays whose RBE was 
generally considered to be higher than that of y-rays. The p-values 
obtained from Welch's t-test between Q^^e ^or each internal 
exposure and the reference condition are also given in Table 1. 
Except for ^^^Cs, the p-values were smaller than the significance 
level; p<0.05. This result indicates that their RBE values are 
greater than 1 in statistically significant. 

RBE for the bystander effect. Our bystander model 
contains four free parameters: a, fi, Mq, and f]. These parameters 
should have a complicated dependence on the biological endpoint 
and the cell type, and the evaluation of their numerical values is 
outside the scope of this paper. Thus, we estimated the RBE for 
the bystander effect for various conditions by arbitrarily changing 
the parameters. 

As examples, the bystander probabilities for exposure to 
electrons from ^^^Cs, a particles from ^^^Pu, and 662 keV photons 
are shown in Fig. 7 for two different parameter settings: 
7Vo= 100,000 and ^ = 0.5 as well as 7Vo= 10,000 and ^ = 0.005. 
These conditions represent the reactive and unreactive bystander 
signals, where approximately 0.001% and 1% of cells, respectively, 
must be initiated in order to fully induce the bystander effect. For 
both conditions, the a parameter was set to 0 to express the 
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Figure 6. RBEs for the Rl-inhomogeneity effect as a function of the cell nucleus radius. The upper, middle, and lower panels are for the 
cases that RIs are localized in cell nucleus, cytoplasm, and extracellular space, respectively. 
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threshold behavior of the triggering probability, whereas the P 
parameter was changed from 1 to 100 Gy~^ to investigate the 
dependence of the RBE on the threshold specific energy. These 
data were similar irrespective of whether the RIs were uniformly 
distributed inside a lattice or had a microscopic localization 
tendency. 



Nearly perfect agreements between the bystander probabilities 
for the exposures to electrons from ^^^Cs and 662 keV photons 
were observed in Fig. 7 for all parameter settings. The 
corresponding data for the exposure to electrons from ^^^Cs and 
^^^I also agreed with the photon data, although not shown in the 
graph. These agreements were attributed to the fact that their 



Table 1. Mean RBEs and their standard deviations (SD) for the track-structure effect estimated from c/ave(y) in combination with 
the Q(y) relationship. 





Rl 


Mean RBE 


SD 


P-value 




1.00 


0.014 


0.40 




1.02 


0.036 


0.028 


131, 


1.02 


0.016 


0.035 




1.49 


0.000 


0.000 




37.0 


0.26 


0.000 



The p-values obtained from Welch's t-test between the average quality factors for each internal exposure and the reference condition are also given. 
doi:1 0.1 371/journal.pone.0099831 .tOOl 
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Figure 7. Bystander probabilities for various irradiation 
conditions as a function of the mean absorbed dose D, Panels 
A and B show the results for the reactive (A/o= 100,000 and ^ = 0.5) and 
unreactive (A/o = 10,000 and ^ = 0.005) bystander signals, respectively. 
doi:l 0.1 371 /journal. pone.0099831.g007 

radiation fields characterized by/i ave(4 were nearly identical to 
one another, as shown in Fig. 4. Thus, the RBE for the bystander 
effect for exposure to ^^^Cs, ^^*Cs, and ^^^I should be very close to 
1 irrespective of the calculation conditions. 

Conversely, the bystander probability for the exposure to cx 
particles from ^^^Pu was significantly higher than that for the 
electron and photon exposures at the same absorbed dose. This 
tendency can be explained by the probability of cell nuclei having 
the specific energy over z in those radiation fields with the mean 
absorbed dose D, f^^Q{>z,D), where their examples are shown in 
Fig. 8. For electron or photon exposures with Z)= 10 |lGy, 
approximately 1% of the cell nuclei were irradiated and their 
maximum specific energy was around 10 mGy. Thus, their 
bystander probabilities were very small, even for the largest 
parameter, i.e., the lowest threshold specific energy. On the other 
hand, only 0.004% of the cell nuclei were irradiated for a 10 jiGy 
exposure to d particles, but most irradiated cell nuclei had a 
specific energy >100 mGy. Thus, the bystander effect is induced 
even by such low-dose exposure for the case of the reactive 
bystander signal. Consequently, the RBE for the bystander effect 
for exposure to ^^^Pu is generally >1, and depends on both the 
model parameters and mean absorbed dose in a complicated 
manner. 

Maximum possible RBE. It is evident from the above 
analyses that the Rl-inhomogeneity effect is the dominant factor in 
determining the RBE for internal exposure to ^^^Cs, ^^*Cs, and 
^^^l. In this subsection, their maximum possible RBE was 
estimated considering the realistic RI localization tendency. 
Figure 9 shows the calculated RBE for the Rl-inhomogeneity 
effect for exposure to electrons emitted from ^"^^Cs, ^"^^'^ 




'Pu 

662keV photon] 



Cell-Nucleus Specific Energy, z (Gy) 



Figure 8. Probabilities of cell nuclei having the specific energy 
over z in various radiation fields, 4ve(>^/^)» Data are for exposure 
to electrons from ^^^Cs, cx particles from ^^^Pu, and 662 keV photons 
with mean absorbed dose D=10)iGy, 1 mGy and 100 mGy, 
respectively. 

doi:10.1371/journal.pone.0099831.g008 

as a function of the RI fraction in cell nuclei. In this calculation, 
the rest of RIs was assumed to be uniformly distributed in the 
cytoplasm and extracellular space. These data were for the largest 
cell size [r^ = 7 jlm, rc = 2 1 jlm), which yielded the highest RBE as 
shown in Fig. 6. The RBE linearly increased with increasing RI 
fraction. 

According to an animal study in Ref [9], 12— 21 % of ^^^Cs were 
localized inside the cell nuclei. If the maximum fraction of ^^'^Cs 
and ^"^^Cs in nuclei of human cells is also 21%, then the maximum 
possible RBE for exposure to electrons from those RI will be —1.1 
and —1.2, respectively. This estimate is quite conservative because 
cell-nucleus mass adopted in this calculation was very small - only 
2% of the total weight. Thus, intranuclear cesium concentration is 
approximately 10 times higher than that in the other structures. 
This value is quite high compared to the results obtained with the 
yeast Saccharomyces [33], where intranuclear cesium concentration 




RI Fraction in Cell Nucleus 

Figure 9. RBEs for the Rl-inhomogeneity effect as a function of 
the RI fraction inside cell nuclei. Data are for exposure to electrons 
emitted from ^^^Cs, ^^"^Cs, and ^^^1, and for the largest cell size 
[r^ = 7 jim, rc = 21 jim). 
doi:10.1371/journal.pone.0099831.g009 
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and cytoplasmic one were nearly equal. In addition, the electron 
contributions to the effective dose were estimated to be only 43% 
and 16% due to the intake of ^^^Cs and ^^*Cs, respectively, on the 
basis of the specific absorbed fractions calculated using the ICRP/ 
ICRU adult male reference phantom [34] in combination with a 
physiologically based biokinetic model for cesium in human body 
[35]. The rest is the contribution from photons whose RBE is 
equal to 1 , and thus, the maximum possible RBE for the intake of 
^^'^Cs and ^^*Cs including the photon contribution were approx- 
imately 1.04 and 1.03, respectively. 

For the exposure to ^^^1, there is no evidence that thyroid cell 
nucleus accumulates iodine, although its subcellular distribution 
has not been thoroughly revealed. In addition, the extracellular 
space such as lumen and blood contains a non-trivial portion of 
iodine [36]. On the other hand, the thyroid hormones containing 
iodine tend to bind to nuclear receptors in cells of extrathyroidal 
tissues such as liver and kidney. However, the organ absorbed dose 
in these tissues due to the intake of ^^^I was much smaller than that 
in thyroid, which accounted for —98% of the effective dose. Thus, 
the assumption that iodine is uniformly distributed inside the cell 
nuclei probably yields an adequate estimate of the RBE due to the 
intake of ^^^I. Even if the intranuclear iodine concentration is 
twice as high as that of the entire system, the RBE would only be 
1.02. 

Conclusions 

The PDs of specific energy inside cell nuclei and dose PDs of the 
lineal energy for a site diameter of 1 |Im for internal exposure to 



^^^Cs, ^^^Cs, and ^^^I, as well as external exposure to 662 keV 
photons, were calculated by performing Monte Carlo particle 
transport simulations with PHITS. The RBEs for the RI- 
inhomogeneity, track-structure, and bystander effects were then 
derived from the calculated PDs. The RBEs for the track-structure 
and bystander effects were very close to 1, owing to the nearly 
identical PDs between the internal exposure and the external 
exposure to 662 keV photons. On the other hand, the RBEs for 
the Rl-inhomogeneity effect largely depended on the intranuclear 
RI concentration and cell size. However, their maximum possible 
RBE was only 1.04 even under conservative assumptions. 
Therefore, it can be concluded from the microdosimetric 
viewpoint that the risk from internal exposure to these RIs should 
be nearly equivalent to that from external exposure to y-rays at the 
same absorbed dose level, as suggested in the current ICRP 
recommendations [22]. 
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